Moduli evolution in the presence of thermal corrections 
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We study the effect of thermal corrections on the evolution of moduli in effective Supergravity 
models. This is motivated by previous results in the literature suggesting that these corrections could 
alter and, even, erase the presence of a minimum in the zero temperature potential, something that 
would have disastrous consequences in these particular models. We show that, in a representative 
sample of flux compactification constructions, this need not be the case, although we find that 
the inclusion of thermal corrections can dramatically decrease the region of initial conditions for 
which the moduli are stabilised. Moreover, the bounds on the reheating temperature coming from 
demanding that the full, finite temperature potential, has a minimum can be considerably relaxed 
given the slow pace at which the evolution proceeds. 



I. INTRODUCTION 



The study of moduli evolution is an active field of re- 
search in the context of the phenomenology and cosmol- 
ogy of string models. At the level of the D — A effective 
theory, which we normally assume to be A/" = 1 Super- 
gravity (SUGRA), moduli are complex scalar fields, some 
of which physically parametrise the size and shape of the 
six or seven original string dimensions that have been 
compactified. It is therefore mandatory that any realis- 
tic model provides, at the end of the day, moduli with a 
non trivial vacuum expectation value at the right scale 
(which is the Planck mass, Mp for conventional small ex- 
tra dimensions), and that the minimum corresponds to 
almost Minkowski space and Supersymmetry broken, in 
order to connect with the Standard Model. Throughout 
the past twenty years there has been steady progress in 
our understanding of the dynamics of moduli, with two 
outstanding problems having to be addressed: the first 
one is the recurrence of a negative (i.e. anti de Sitter, 
AdS) vacuum energy for every model in which moduli 
were successfully stabilised, while the second one is the 
fact that these potentials are so steep along certain direc- 
tions that, from the dynamical point of view, it looked 
impossible that given any initial conditions away from 
their minimum, the moduli would end up in it. This is- 
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sue was first pointed out by Brustein and Steinhardt [l[ 
and is commonly known as the problem of the " runaway 
dilaton" . 

Concerning the first problem, namely the recurrence of 
AdS solutions within all successful attempts at stabilising 
moduli, recent developments in the context of flux com- 
pactifications in type IIB string theory 0] have opened up 
new ways of trying to achieve either Minkowski or de Sit- 
ter (dS) vacua. Iiiparticular, the mechanism presented 
by Kachru et al. Q (KKLT from now on) realised this 
by adding D-terms to a SUSY-preserving, AdS F-term 
vacuum. Although not entirely correct in the context 
of Supergravity (see 0, Q for some criticism, and 0| 
for proposed solutions), its main features have triggered 
an enormous amount of interesting work, and subsequent 
progress in this field through the past years (in particu- 
lar, for explicit string realisations, see [3, Q). 

Also on the topic of runaway moduli, in general, 
substantial developments have taken place in the past 
decade. For a range of initial conditions the problem can 
be alleviated considerably by considering a background 
perfect fluid which decelerates the fields and prevents 
them from passing the barrier dividing the physical vac- 
uum from the runaway one [l^, [ll|, 0, [Hj [3 ■ 
Thermal corrections, however, are potentially danger- 
ous as they modify the shape of the potential and, 
at high temperatures, the physical vacuum is entirely 
lost [itI . [Tsj . However, there is a potential limit to 
the extent to which this argument can be used. In the 
very early Universe as the value of the Hubble param- 
eter, becomes close to the Planck scale, scattering pro- 
cesses are unable to establish thermal equilibrium be- 
cause they do not have sufficient time compared to the 
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expansion rate of the Universe [l9|. In this era thermal 
corrections arising out of these scatterings would not be 
present, and so the physical vacuum would not be desta- 
bilised. In the context of the Standard Model Enqvist 
and Sirkka considered the thermalisation of a hot QCD 
gas in the early universe, and calculated the critical tem- 
perature above which the Universe can not thermalise 
to be Tcrit = 3 X 10^^ GeV HO]. In Ref. [2l| the authors 
used this argument in considering racetrack inflation and 
assisted moduli stabilisation. 

In this work we present a detailed quantitative anal- 
ysis of the effect of these thermal corrections on the re- 
gion of initial conditions leading to stabilisation of the 
moduli, for a given SUGRA model. Given that we are 
working in regimes beyond the standard model, where 
we do not have a proper handle on the conditions un- 
der which we will move out of thermal equilibrium, we 
adopt two approaches. The more theoretical is to accept 
that the Universe may have been in a period of thermal 
equilibrium close to the Planck scale and investigate the 
impact of thermal corrections in those regimes. The sec- 
ond is to take seriously the QCD bound, and investigate 
the impact of the corrections in these lower temperature 
regimes. In both cases we will conclude that thermal cor- 
rections do decrease the area of the region of stabilisation 
by a similar relative amount. This effect from the ther- 
mal coupling will act on top of the previous results, where 
the absolute size of the stabilization region decreases for 
a smaller initial density of the background fluid pj"'* [l^ . 



II. EQUATIONS OF MOTION 



We will be treating Oq and 02 as variables which can be 
varied in order to test the results for a range of possible 
values of Nc and Nf. We can obtain the energy density 
and pressure of this thermal fluid as p,. = —F and p,, = 
—Pr + Tdpr/dT, hence 



(6) 



where r = a^ja^. 

In principle, a second component of relativistic parti- 
cles that only interacts with the moduli fields gravita- 
tionally can also be present in the dynamics. Assuming 
both components of radiation to be in thermal equilib- 
rium, the non interacting radiation will have an energy 
density given by pB = 7r^(7*T'*/30 where is the num- 
ber of effective degrees of freedom at the temperature T 
of the thermal fluid p,.- In this case, the equations of 
motion above arc still valid with ao being replaced by 
ao ^ flQ — 7r-^(7*/90. The effective |ao| can then be very 
large , so that r w 0, effectively washing out the effects 
of the thermal corrections on the moduli evolution. 

The effective scalar potential for the moduli in four 
dimensional A/" = 1 Supergravity is given by 



V = e^(K^'^DiWDjW - 3IUIU) 



(7) 



where if*-' is the inverse Kahler metric and DiW = diW+ 
diKW is the Kahler covariant derivative acting on the 
superpotential. In general, the Kahler potential AT is a 
function of the real parts of the fields and, for most string 
compactifications, acquires the form 



A' = -^ln($, + $,) 



(8) 



In this work, we will be studying two string theory 
models that can be described by a four dimensional N = 
1 effective Supergravity theory with action of the form 



S 



where 



-R-C^ + F{g,T)] , (1) 



U^-K.jd^^'d^'^ -V , (2) 

and Kij = d'^K/d^^d^^ is the Kahler metric; are 
complex moduli scalar fields; V{^) is the scalar potential 
and G is the 4-dimensional Newton constant. The free 
energy F[g,T) acts as a Lagrangian density of matter 
fields. For a SU(A^c) gauge theory with Nf multiplets at 
high temperature T, the free energy has a perturbative 
expansion in terms of the gauge coupling g = 
where $h denotes the real part of the moduli fields <f>. 



F{g,T)^{ao+a2g^)T\ 
and the parameters ao and 02 are given by 



ao 
02 



24 



{n! + 2N,Nf - 1) , 



64 



{N^-l)iN, + 3Nf). 



(3) 

(4) 
(5) 



where the sum is understood over all moduli <f>i. As for 
the superpotential, which encodes the dynamics of these 
fields, a combination of flux terms (normally polynomials 
in the different moduli) and non perturbative effects (in- 
stantons and gaugino condensation being the most well- 
known ones [23, |2^, [U, i25i|) will provide the potential 
with a non trivial vacuum structure. For the purposes of 
showing the effects of thermal corrections, we use the toy 
models of KKLT Q and Kallosh and Linde [H (KL from 
now on) , given that they capture the essential features of 
the cosmological problems usually attributed to moduli. 

The equations of motion follow from the variation of 
the action ([T]). Considering homogeneous fields evolving 
in a spatially flat Friedmann-Robertson- Walker space- 
time background, the equations of motion for the com- 
plex fields yield 



rpr 



3(l + rg2) 



K^'djg' 



(9) 



where = d^^dt, djV = dV/d^^, and the connection 
on the Kahler manifold has the form 



(10) 
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In addition, the Hubble rate H = a/a, where a{t) is the 
scale factor of the Universe, is given by the Friedmann 
equation 



Mp\p^+Pr)^Mp\K,j^'^ +V + Pr) , (11) 



with Mp 2 = 8ttG (or Mp = 2 x lO^^ GcV) and ^ 
• -J 

Kij^^ ^ +V and pr are the energy densities of the evolv- 
ing moduli fields and of the thermal fluid, respectively. 
In what follows we set Mp = 1. 

We need to understand now how the temperature T 
relates to the values of the fields and the scale factor of 
the universe. To this end we note that the equations of 
motion for the scalar fields can be rewritten as 

1 -1 -J 

-SHip^ + p$) + -^^pr^' + -djPr'^ , (12) 

where the pressure of the moduli fields is defined as ~ 
■ ■ -J 

Kij^^ ^ — V. By requiring conservation of the total 
energy density p ~ p<s, + Pr we must have 

1 -1 

Pr = -AHpr - -d,pr<^' - -djPr<^ , (13) 

which, upon integration, gives a solution for pr of the 
form 

Comparing with Eq. ^ , the evolution of the temperature 
can be seen to be 



T = Ti,- 



l + r(72($^) 



(15) 



It is worth splitting the equations of motion for the 
complex scalar fields into those for their real and imagi- 
nary parts 



$^ -t- 3i/$^, + r;.,($j,$5^^ - $^,$J) + ^-K^'d.^v (16) 



1 
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^]+3H<i>]+T),,m<P^i^+^'j^^'j) + -K'^dj,V = , (17) 

where now ($}) refers to the real (imaginary) part 
of the scalar fields and {djj ) are used to denote the 
derivative of the potential with respect to the real (imagi- 
nary) parts of the fields respectively. Wc note that, given 
an initial value of the energy density of the thermal fluid, 
p^"'*, the equations of motion and the Friedmann equa- 
tion only depend on the ratio r = a2/aQ but not on the 
specific values of ao and 02. 

Given Eqs. ^ and wc can define an effective 
scalar potential for the fields. 



Feff($) = F($) 



1 rprg^ 
3 1 -l-ro 



(18) 



up to a constant. With a coupling constant of the form 



c 



(19) 



where c is a constant, it is clear that, at high temperature 
(large pr), the effective potential can look very different 
from V and, in particular, it can be devoid of a mini- 
mum, limiting the stabilisation of the moduli or the max- 
imum t emp erature allowed. This is the problem raised 
in Ref. 



III. KKLT MODEL 

The possibility of finding de Sitter vacua in string the- 
ory with a stabilized volume modulus, a, was put forward 
in Ref. and has been widely adopted in subsequent 
work. The key ingredient was to consider the combina- 
tion of non perturbative effects and an additional flux 
term in the superpotential 

W = Wa+ Ae-""" , (20) 

which, together with the usual Kahler potential 



K = -31n(cr + a) 



(21) 



defines the F-term of the SUGRA potential, see Eq. ([7]). 
It has been known for many years now that, in this con- 
text, it is possible to stabilize cr, although giving rise to 
an AdS vacuum. As pointed out in Ref. [3|, however, if 
contributions from either anti-D3 or D7 branes are in- 
cluded, an additional D-type term of the form 



(22) 



is generated, where we write cr = an + iaj. By suitably 
tuning the value of C one can move to a de Sitter - or 
even Minkowski ~ vacuum. The scalar potential for a 
has an extremum in ai for aaj = tt-tt, with n an integer. 
Depending on the sign of Wq cos{aai) this can be either 
a maximum or a minimum. 

In this work we are interested in studying the cosmo- 
logical evolution of the field a as it rolls towards its min- 
imum. Previous ai ialy sis addressing the same issue were 
published in Refs [li,[il|, however, without taking into 
account the effect of the thermal corrections. Here, we 
compare the profile and area of the region of initial po- 
sitions of the field, CTinit that leads to stabilisation when 
the field rolls in the presence of a background perfect 
fluid with the same region of stabilisation when the field 
evolves in the presence of thermal corrections. We illus- 
trate that comparison for the KKLT model in Fig. [T] and 
Fig. [2] for two values of the initial energy density of the 
thermal fluid, pj,"'*. The first corresponds to the very 
high energy regime p""* = 10~'*Mp, where the second 
corresponds to a value which satisfies the condition for 
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thermal equilibrium obtained in [20] for |ao| < 0(100), 
namely pj"'' = In both cases the stabilisation 

region is shown with thermal corrections (shaded areas) 
for given values of the ratio r = 02/00 against the sta- 
bilisation region with just a perfect fluid of radiation (or 
equivalently, with r = 0) for the same initial pr (solid 
black line). 

It proves convenient to work with the canonically nor- 
malized field (j) = ^3/2 In cr/^, instead of an itself. In the 
lower right panel of the figures we draw the ratio between 
the area of the two stabilisation regions — i.e. the one 
with thermal corrections (Ath) and the one without {A) 
— against \r\. 

We note that, because the energy density of the ther- 
mal fluid must be positive definite, and taking the effec- 
tive 4D Yang-Mills coupling on the D7 brane = 47r/(Tij, 
we must ensure that 



r = -0.7 



r = -1.3 



1 

47r 



/2730 < ^ < 



(23) 



This means that we need to impose a larger lower bound 
on the initial values of </> as we increase \r\. In practice 
this bound is only effective for large initial p™'* since the 
stabilisation region decreases naturally for lower initial 
temperatures. 
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FIG. 1: The panels with the shaded regions show the regions 
of initial conditions (0, 01) - in Planck units - that lead to 
stabilisation of a at the minimum of the potential for the 
KKLT model in the presence of thermal corrections. The 
solid black line corresponds to the region of stabilisation in 
the presence of a perfect fluid (r=0) with the same initial 
energy density p!-"'*. The lower right corner shows the ratio 
of the areas of these regions against the ratio r — a2/ao- We 
have set pj"'* = 10'^ M^. 



Returning to Figs. [THll we show the regions of ini- 
tial conditions that lead to the field stabilising in the 
minimum of its potential. We used the same values of 
the parameters for the model as in Ref. [l^ . namely 




FIG. 2: Same as Fig.[T]but setting p^"" = 10""A/p. 



A = 1.0, a = 0.1, C = 3 X 10^26^ ^^^^ negative 
(with cos{aai) — 1) such that the minimum at cr/ = 
is supersymmetric. Comparing the two figures, we can 
see that the stabilisation region for the lower initial tem- 
perature is considerably smaller, a result that was shown 
previously for the evolution without a thermal coupling. 
The effect from the thermal coupling can be more clearly 
seen in the lower right plots, where the ratio of the two 
stabilisation areas is shown against \r\. We see that in 
both cases increasing the strength of the thermal coupling 
(that is, the relative size of 02) can effectively eliminate 
the stabilisation region. On the other hand, the curves 
for the two different choices of initial p™'* are very similar, 
meaning that the relative effects of the thermal coupling 
seem to be independent of the initial value of the back- 
ground. In both cases, a value of \r\ ^ 1 implies that the 
stabilisation region is decreased by less than 50%. 

Using Eq. © with oq = —100, we can see that our 
choices for the inital background energy density cor- 
respond to initial temperatures of T ~ 10~^Mp and 
T ~ 10^^ Mp. However the dilaton potential only de- 
velops a minimum at a temperature Tcrit ~ 10~^Mp, so 
that we can have stabilisation of the dilaton even when 
the evolution starts at a temperature where the minimum 
does not exist. 

This behaviour can be easily understood. Thermal cor- 
rections act as a background fluid in that they bring an 
extra contribution into the Friedmann equation. Though 
at high temperatures the structure of the effective poten- 
tial is spoiled (that is, we have no minimum), the extra 
contribution reinforces the frictional term in the equation 
of motion and forces the field to slow down its evolution 
as when in the presence of a non interacting perfect fluid. 
Therefore, for suitable initial conditions, the field can ap- 
proach the value corresponding to the minimum of the 
T ~ potential, {4>)t=o, when the thermal corrections 
have already become negligible, i.e. the temperature has 
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decreased below the critical temperature for which the 
minimum is created. In Fig. [31 we show, in the left panel, 
the evolution of and pr as a function of the tempera- 
ture for the model used in Fig. [U having singled out five 
values of the temperature. In the right panel we can see 
the profile of the effective potential for these five values 
of the temperature with the corresponding value of the 
field at that particular time. We can sec, in position 3, 
at T = 10~^Mp, that the field is evolving in an effective 
potential without a minimum but that, in position 5, the 
effective potential has recovered the minimum and the 
field has been trapped in it. 

We observe that, initially, the large kinetic energy 
forces the field to evolve to the flatter region of the poten- 
tial where the field effectively freezes due to the large fric- 
tional term in the equation of motion given that p,. ^ 
holds. As the universe expands the thermal corrections 
continue to decrease in magnitude and the field reenters 
the steep slope of the potential which admits a scaling- 
like regime of evolution i.e., the field restarts rolling with 
a kinetic energy that is proportional to the potential en- 
ergy. It is this property that prevents the field from gain- 
ing kinetic energy and going over the barrier that sepa- 
rates the physical minimum from the runaway one. In 
other words this is why the modulus does not get driven 
to larger values more quickly. For an analytic descrip- 
tion of this stabilisation mechanism due to scaling, see 
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FIG. 3: In the left panel we show the evolution of the en- 
ergy densities of the field, (solid line) and of the thermal 
fluid, pr (dashed line) in the KKLT model, as a function of 
the temperature T (all in Planck units) . In the right panel 
we show the profile of the effective potential at five differ- 
ent values of the temperature (solid lines) and the position 
of the field at these different stages (circles), as a function of 
(f) (also in Planck units). We have set r — —0.2, ao — —100, 



r = -0.4 



r = -0.7 






FIG. 4: Same as in Fig. [T]for the Kallosh-Linde model with 
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IV. KALLOSH-LINDE MODEL 

The Kallosh-Linde model (KL) [2^ generalizes the 
original version of the KKLT model by admiting two 
components in the superpotcntial 



W = Wo + Ae~ 



Be 



(24) 



10" 



3.2 and ai = 0. 



A particular example was investigated in a previous pub- 
lication The parameters of the model were set to 
A= 1, B = -1.5, C = 0, a = 27r/f00, /? = 27r/99 and 
Wo such that there is a supersymmetric minimum with 
zero cosmological constant, i.e. Wq is such that both W 
and Fa = K^W + Wa vanish at some cr_R = (Tcrit, 0-7 = 0. 
Furthermore, there are a series of supersymmetric, AdS 
minima. In Figs. |4][5] we show how the stabilisation re- 
gion with thermal corrections varies with r and, for com- 
parison, we plot it against the stabilisation region when 
a perfect fluid of radiation is present (or equivalently, 
when r = 0). As for the KKLT model, we see that for 
both energy regimes the profile of the stabilisation re- 
gions is practically unmodified for small values of |r| and 
that their area decreases as \r\ becomes larger, eventually 
vanishing for sufficiently high values of \r\. Furthermore, 
we see again that this effect of the thermal coupling on 
the stabilisation region is similar for the two initial values 
of pj,"'* in relative terms, though a smaller initial value of 
the initial energy density in the background makes the 
stabilisation areas smaller in absolute terms. The ther- 
mal coupling is slightly more effective in reducing the 
stabilisation area in this KL model. 



V. DISCUSSION 

We have seen in the two models investigated that the 
inclusion of thermal corrections in the moduli fields evo- 
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r = -0.4 ,10-= r = -0.7 




FIG. 5: Same as in Fig. [T] for the Kallosh-Linde model with 

lution can have a very similar effect to a perfect fluid 
of non interacting radiation, in the sense that their con- 
tribution increases the friction, decelerating further the 
fields thus helping in their stabilisation at the desired 
minimum. We saw that the relative effect of the thermal 
coupling in the stabilisation area of the dilaton is nearly 
independent of the actual initial value of the background 
energy density. As we increase the value of |r| (that is, 
we increase the relative value of 02), the stabilisation re- 
gion becomes smaller, and eventually disappears. Note 
however, that in both models, values of |r| ^ 1 will only 
decrease the area for at most 50%. For reference, to 
have |r| > 1 we need Nc > 17, whereas for \r\ > 1.5 
we need Nc > 26. The values of \r\ obtained strictly 
from the thermally coupled fields only gives us a lower 
bound on the stabilisation area. In principle, we will also 
have fields in the background without a thermal coupling 
to the dilaton in the background. In this case, wc will 
have an effective value for |7'| <^ 1, reducing consider- 
ably the global effect of the thermal coupling. Namely 
for \r\ < 10~^, the difference between the two scenarios 
is of less than 1%. On the other hand, it should also be 
clear that a larger value of the background energy den- 
sity leads to a wider region of stabilisation, as the field 
can enter a scaling regime earlier. 

The existence of a minimum in the potential is usually 
seen as a necessary condition for the stabilisation of a 
moduli field, from which upper bounds to the reheating 



temperature can be obtained. We have seen that, if we 
allow for the field to evolve, the stabilisation is also de- 
pendent on its initial condition; and a minimum in the 
potential need not be present initially. From Eq. ([5]) we 
get that for \r\ ^ 1 the relation between the initial value 
of the temperature and energy density is approximately 

which means that, for values, Oq « —100 and pj"'' ~ 
lO-^M^ or p™* w 10-"M^, the initial value of the 
temperature can be as large as 10~^A/p or 10~'*A/p, re- 
spectively. These are above the usual upper lim it g iven 
for stabilisation, T < 10~*Mp suggested in Ref. [Igl ob- 
tained by requiring the minimum to appear in the effec- 
tive dilaton potential. 

Even though we do consider sectors different from the 
(MS)SM as the source for the thermal bath, the simplest 
scenario would be to assume that the SM fields would 
eventually enter thermal equilibrium at the same tem- 
perature, assuming a smooth evolution of all the back- 
ground fields. Of course this does not need to be the case, 
and in that instance a different later background with a 
higher temperature that interacts with the dilaton could, 
in principle, destabilize the field again. Assuming the 
field to be already in its minimum, this would have to be 
at a sufficiently high temperature to effectively remove 
the minimum. 

Though the mechanism here described seems very gen- 
eral, it would be interesting to evaluate its robusteness 
for other scenarios such as the KKLT model coupled to 
a Polonyi field studied in Ref. [l^. One other aspect 
is that to large temperatures correspond large thermal 
fluctuations which may give rise to large spatial inhomo- 
gcneitics in the moduli. A quantitative investigation of 
these effects and their cosmological implications deserves 
a complete study. 
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